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This paper deals with the concept of adiabaticity for fully quantum mechanically cavity QED 
models. The physically interesting cases of Gaussian and standing wave shapes of the cavity mode 
are considered. An analytical approximate measure for adiabaticity is given and compared with 
numerical wave packet simulations. Good agreement is obtained where the approximations are ex- 
pected to be vàlid. Usually for cavity QED systems, the large atom-field detuning case is considered 
as the adiabatic limit. We, however, show that adiabaticity is also vàlid, for the Gaussian mode 
shape, in the opposite limit. Effective semiclassical time dependent models, which do not take into 
account the shape of the wave packet, are derived. Corrections to such an effective theory, which 
are purely quantum mechanical, are discussed. It is shown that many of the results presented can 
be applied to time dependent two-level systems. 
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I. INTRODUCTION 

The concept of adiabaticity dates back to the early 
days of quantum mechanics and has, in spite of its 
simplicity, drawn attention ever since. In this pa- 
per we investigate adiabaticty within cavity quan- 
tum electrodynamics (QED) models, namely for ex- 
tended Jaynes-Cummings models 0. A two-level 
atom interacts with a quantized cavity mode, whcrc 
both atòmic center-of-mass momentum and posi- 
tion are taken into account. The adiabatic theorem 
states that, in the limit of an infinitely slowly chang- 
ing Hamiltonian H (t) time t may be seen as a param- 
eter and H(t) can be instantaneously diagonalizable 
and no population transfer will take place between 
the instantaneous eigenstates [2| . If the change gov- 
erns a closed loop in 'parameter' space, it has been 
shown that the eigenstates do not only accumulate 
the Standard dynamical phase, but also a geomètric 
one . Adiabatic methods have also been suggested 
for realizing quantum information processing, caus- 
ing the parameters to be stable against fluctuations 
U- Recently, an inconsistency of the criteria for 
the theorem has been proposed || , which has led to 
new insight into the theorem and several comment- 
ing papers |(| . The adiabatic theorem is usually as- 
sociated with an explicitly time dependent problem. 
However, the phenomena has a broader spectrum 
of application, for example: adiabatic elimination 
of infernal states 0, II] and the Born-Oppenheimer 
approximation 0, and within wave packet propa- 
gation. Adiabatic elimination gives effective models 
when the variations of some slowly varying quanti- 
ties are neglcctcd from the equations of motion. For 
wave packet propagation, which will be considered 
in this paper, the rate of change is governed by the 
motion of the wave packet. This is easily seen for 



particle with mass m and kinetic energy much larger 
than the interaction energy; it is often legitimate to 
replace momentum and position in the Hamiltonian 
by their classical counterparts, p — ► mv and x — > vt. 
In this way an effective time dependent model is ob- 
tained |lCt which does, however, not take into 
account the shape of the wave packet. It is known 
that, for small kinetic energies, treating x and p as 
classical variables may not be justified, and a full 
quantum analysis must be carried out |12L Il3| . 

We introduce a dynamic and adiabatic frame in 
sections lïïl and lïïïl and give a criterion for adiabatic 
evolution of the wave packets in section IIVI This 
measure is approximated to first order and checked 
with numerical wave packet propagation. Good 
agreement is obtained, when expected according to 
the approximations. Higher order corrections, for 
example when the shape of the wave packet is taken 
into account, are discussed in section IVÏÏ1 The gen- 
erally associated adiabatic limit |A| ^> G(x), where 
G(x) = ^/ng{x) is the atom-ficld coupling constant 
multiplied by the square root of the photon number, 
reproduces adiabaticity in our model. It is, however, 
found that the opposite limit, |A| <C G(x), can also 
give an adiabatic process. This is counterintuitive 
from what is formally believed, and is therefor dis- 
cussed in detail in the paper. Two special cases of a 
Gaussian and a standing wave mode shape are con- 
sidered separately in subsections IV Al and IV Bl for 
a theoretical analysis and in IVI Al and IVI Bl for nu- 
merical simulations. In the former Gaussian model, 
with |A| <g; G(x), an adiabatic evolution is obtained 
for all x, while for the standing wave mode with 
wave length A, the adiabaticity breaks down around 
the points x = tm/q, where n = —2, —1, 0, 1, 2, ... 
and q = 2tt/X. As the approximated adiabaticity 
criteria have identical form as the adiabaticity cri- 
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teria obtained in time dependent models, much of 
the analysis in this paper also applies to such sys- 
tems. This is discussed in section IVIIII where we 
also give a systematic way to derive a semiclassical 
time dependent Hamiltonian from time independent 
fully quantum mechanical ones. 



II. THE DYNAMIC FRAME 

We can easily write down the Hamiltonian de- 
scribing a single two-level atom passing an electro- 
magnètic field confined to a cavity. In addition to 
the customary Jaynes-Cummings terms, we have to 
add the energy of the free atom and the spatial vari- 
ation it feels from the caity. With Tí = 1, we obtain 

v 2 £1 
Hi = ^— +wa T a+ -0-3 + H 2 (x), (1) 



we find with 



(2) 



where the cavity mode interaction is 

= g(x) {a<j + + a'cr~) . 

Here Í2 is the energy separation of the two-level sys- 
tem, ui is the energy of the photon created by a\ 
and the a : s are the Pauli matrices. The excitation 
number operator 



N = a)a+ -<t 3 



(3) 



commutes with the Hamiltonian and thus we find 
separation of the blocks of the state space labelled 
by the excitation number 

3»(p) = 3+(p) | +,n- 1) + S_(p) | -,n). (4) 

The translational degree of freedom is found in the 
coefficients $±. 

We may, however, separate an arbitrary múltiple 
A of N because of its commutation with the Hamil- 
tonian. We dcfine the states 



$" = cxp(-iAiVí)$ ri 
cvolving with the Hamiltonian 
H = H 1 - KN 



(5) 



(G) 



= p 

2m 



+ {uj-A)aTa+ ^a 3 + V(x) 



These are the generalisations of going to a rotating 
frame. We now have two important cases: 

Case 1: A = uj. 
Case 2: A = fi. 

If we write the equations for the state in the 
general form 




'dl 









+ 


~e+ G ' 








2m 






G e- 







(7) 



the results 



and 



A = D,-lj ; G = g(x)y/n 



Case 1 : 



Case 2 : 
e + =-A(n- 1) 
e_ = —An. 



(8) 



(9) 



(10) 



III. THE ADIABATIC FRAME 

If we now define the two-level potential from (Q) 



V(x) 



e+ G 
G £_ 



we may introduce the unitary transformation 
U 



cos u sm t 
— sin 9 cos ( 



such that the adiabatic potential 

uvu ] = w 



A+ 
A_ 



is diagonal. However, since 

U(x)pU f (x) =p-a 2 d6, 

where 

_ d6(x) 



89: 



dx 



(11) 



(12) 



(13) 



(14) 



(15) 



the transformed Hamiltonian is non-diagonal. Ex- 
plicitly we have 



H — H a d + H cor , 



with the adiabatic part 

,,2 



H, 



ad 



p 
2m 



A + 
A_ 



(16) 



(17) 



and the non-adiabatic correction part 
1 



2m 



{dOf -2{d9)p-d 2 9 
2(89)p + d 2 9 (09) 2 



(18) 

The transformed Hamiltonian is given in the adia- 
batic basis 



| t)n = cos9\n, +) +sin%+ 1,-), 
| í) n = -sin 9\n,+) + cos6\n + 1,-), 



(19) 
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and for a given n, a general state is written, in the 
adiabatic basis, as 

9(x,t)=* r (x,t)\ T)„ + %(M)I |) n . (20) 
The eigenvalues are 



Ah 




G 2 , (21) 



where 

Ae = (22) 
and the diagonalisation gives the result 

2G 



tan 26» = 



Ae : 



(23) 



The interesting observation is that in both Case 
1 and Case 2, we find 



Ae 



(24) 



which shows that the energy splitting l|21(l in the 
adiabatic frame becomes the same. On the other 
hand, the splitting is centered on different positions 
because 



= 0; 



(Case 1) 



= -A(n- i) ; (Case 2). 
For large detuning A, we find from (|21|l 



A H 



A G 2 

£ i — i — . 

2 A 



(26) 



As expected, the levels are pushed apart for both 
signs of the detuning A. In this limit, the atòmic mo- 
tion is taking place on the effective potential surface 



V eff = —^~ n - 



(27) 



We see that the spatial dependence of the mode 
structure is enhanced by the number of photons 
in the cavity mode. This limit is the one usu- 
ally intended in the discussion of adiabatic following 
through the cavity-induced potential. 

When the detuning goes to zero, A — > 0, the adi- 
abatic potentials follow the shape of the mode func- 
tions, (|21|) becomes 



A± = £++ c . £ ± 9(x)Vn. 



(28) 



However, formally the adiabatic diagonalisation 
can be performed for all vàlues of A, and we want 
to enquire whether the result (|21|) has any validity 
beyond the customary large detuning limit? 



IV. VALIDITY OF ADIABATICITY 

The correction introduced by the adiabatic ap- 
proximation follows from the transformation of the 
kinetic energy due to the position-dependent diag- 
onalizing operator, seen in eqs. I|14|l and l|15|l . We 
observe that if the angle 9 were equipped with in- 
dependent dynamics, this would be analogous to a 
gauge field situation. The adiabatic approximation 
holds when H ac i dominates the evolution compared 
to H cor , thus when the distance between the adia- 
batic eigenvalues is much larger than the off-diagonal 
term. This gives a condition for adiabaticity accord- 
ingly 

|(A + )< - (A_)i| » ^ \(2(d0)Ph + (9 2 0)i\ ■ (29) 

Here we have averaged over the adiabatic wave 
packet \E r " d (x, í) defrned as the initial wave packet 
propagated by the adiabatic Hamiltonian l|17|) 



*F(x,t) = exp 



P 

2m 



$>i(x,0), i=i, t, 
(30) 

where it is understood that A + corresponds to the 
i =| and vice versa for A_. Hcnce for any operator 
h(x,p): 



(25) (h(x,p)) i =ldx*f*(x,t)h(x,p)*f(x,t), i=l,h 

(31) 

For a more general initial state, any linear combina- 
tion of | |) and | f) states (n — 1), we define the 
adiabaticity parameter as 



1 



(2(d6) P ) l + (d 2 9) l \ , ATO \(2(d0)ph + (d 2 9) 



■N' 



|(A + )|-(A_)i| i |(A + ) T -<A_) T| y 

where N® is the initial population in the state \i), 
i =1, |. Thus At <C 1 implies an adiabatic evo- 
lution for the particular initial state. In time de- 
pendent systems, the adiabaticity constrain is ob- 
tained directly by instantaneous diagonalization of 
the Hamiltonian. In our model, however, this is not 
the case; the diagonalization gives ^ ad (x, 0)j, but its 
time evolution is obtained from l|30|l . 

To obtain some insight into the adiabaticity pa- 
rameter, some approximations are in order. It is 
reasonable to assume that the main non-adiabatic 
contributions, due to the evaluation with the Hamil- 
tonian H cor (I18|) . comes from the off-diagonal term 
d6p. We will therefor neglect the term d 2 9 at the 
moment and discuss its influence in section IYlII For 
a particle moving with a high velocity, we make the 
assumption that p « po for the two wave packets 
^2 \(x), where po is the initial velocity. The posi- 
tion is thus given by x = pot/m+xo (xq is the initial 
position and is taken to be the same for both states). 
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We return to this approximation in section lVHI To 
evalute (|15|l . we take the derivative of i|23|) to obtain 



296» 
cos 2 29 

This gives 



= (1 + tan 2 26») 289 = 



2dG 



39 = 



AÒG 



A^/ndg 



A 2 + 4G 2 A 2 + Ang 2 ' 



(33) 



(34) 



Summarizing the above and using eq. (|32|1 . we find 
the approximate expression 



A°(x) = 



Po\ A^/ndg 
(A 2 + 4n<? 2 ) 3/2 



(35) 



As it stands, A°(x) is given as a function of the pa- 
rameter x, which is, however, equivalent to express- 
ing it as a function of time t by the substitution 
x = pot/m + xo. The expression has got two 
limiting behaviours: For | A |3> 2g^/n we obtain 



A°(x) 



Po\ s/ndg 
m 



(36) 



This gives the condition for the validity of the or- 
dinary adiabatic result l|27|) . In the opposite limit 
| A |< 2gy/ïï, we find 



A°(x) 



'Po\ Ay/ndg 
m ) 8ng 3 



(37) 



This result is remarkable because it favours a high 
level of cavity excitation and it disappears for van- 
ishing dctuning. This sccms to indicate that the 
adiabatic approximation may also be vàlid for reso- 
nance between the àtoms and the cavity mode. How- 
ever, in order to see how this comes about, we need 
to investigate the behaviour in some special cases. 



V. SPECIAL CASES 

A. Gaussian mode shape 

If the àtoms traverse the mode structure in a 
transverse direction of a Fabry-Perot cavity, the am- 
plitude seen is of Gaussian type. Thus we set 



.4 



2ira 



exp 



x 

2^ 



(38) 



For this coupling, the adiabaticity criterion A°(x) <C 
1 reads 



A^a. 



(-)- 

\mJ ( A 2 +4ng 2 (x)f /2 



(:)«« 



< 1. 



(39) 



For large detunings, because x < a, this becomes 
approximately 



(-) 



aA 2 



< 1. 



(40) 



Like (|36(l this is the ordinary condition for adiabatic 
motion; in particular, it is less vàlid for high level of 
excitation. 

In the opposite limit, A —> 0, we find 



1 » 



(-) 



A 



8ng 2 (x) 



Po 



exp 



(41) 

This seems to vanish for A = but, on the other 
hand, it blows up for x > a. 

In deriving (|4l·|l we assumed |A| <C ^Jng{x), which 
will clearly break down for vàlues on x large enough, 
say x > a. This is exactly when the adiabaticity pa- 
rameter A°(x) seems to blow up according to (|41|) . 
For such vàlues this approximation is not likely to 
be vàlid, and we need to explore the bahaviour of 
A°(x) for large |a;| and small A in more detail. Is 
A°(x) approaching smoothly or discontinuously 
when A — > 0? In fig. we plot ^4°(a;) as func- 
tion of x and A (starting from A = 0.0001). In the 
following we consider Case 1: A = u>, and use sealed 
parameters such that H = m = 1 and put the photon 
number n — 1 in all of the numerical calculations. 




- 300 - 200 



FIG. 1: The adiabaticity parameter A°(x) as function 
of x and A in sealed dimensionless units (starting from 
A = 0.0001) for a Gaussian coupling 138H with amplitude 
A = 1, the pulse width a — 50 and the initial momentum 
Po = 10. 

We note that ^4°(a:) has two màxima for a given 
detuning A, symmetrically situated around x = 0. 
The value of this maximum is clearly increasing for 
small detunings, but it is also moving towards |a;| 3> 
a. Thus, in the very limit of vanishing detuning, 
the adiabaticity parameter approaches zero for all 
x. This is more clearly seen in fig. [21 which shows 
the location of the maximum (x > 0) as function of 
A for three different couplings A. 
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FIG. 2: The location of the maximum max{.4 (a;)} as 
function of the dimensionless parameter A for three dif- 
ferent couplings. The pulse width is again a — 50 and 
the initial momentum po = 10. 



Note that the màxima are located at x = ±o for 
large detunings, which might be expected since that 
is where |dg(a;)| is maximal. For small detunings this 
is not the case, and a numerical check shows that 
for A < y/nA/10 the màxima deviate from x = ±o. 
Thus, for small A and within these approximations, 
this seems to confirm that a large coupling or a high 
excitation n gives an adiabatic behaviour, which is 
contrary to the opposite case of A large. 



B. Standing wave mode shape 

If the cavity field contains many spatial periods, 
the atom will sec the potential modulated by several 

periods A = \ ■ The mode coupling then becomes 

g(x) = As'mqx, (42) 
which gives the condition for adiabaticity 



A°(x) 



qAy/ïïA cos qx 



(A 2 + AnA 2 sin 2 qx) 



3/2 



< 1. (43) 



In the ordinary adiabatic limit of large detuning, this 
becomes 



'Po\ qy/nA 
,m) A 2 



< 1. 



(44) 



We find, as we expect, the analogue of 1)41) |): the 
result improves with increased detuning and lower 
cavity excitation. 

The opposite limit A — * is more interesting and 
more subtle to analyze. Cavity QED systems oper- 
ate in two different regimes; the optical fl4| and mi- 
crowave regime [Ï5|. In the former case, the atòmic 



wave packet may extend over several periods of the 
mode shape, and therefore sees a smeared periòdic 
potential, whilc in the latter case A is assumed to be 
larger than the width A^ of the atòmic wave packet. 
This results in very different adiabatic evolution in 
the small detuning limit, which will be seen in the 
next section. 

The adiabaticity parameter A°(x) is a periòdic 
function of x and we may study it around one of the 
màxima, say x = 0. For x ~ we find 



A(x) 



Po 



\C\ 



2mA (i + C 2 x 2 ) 3/2 ' 



(45) 



where C — Iq^fnAj A. Note that this also gives 
the approximated adiabaticity parameter A°(x) for 
a linear coupling g(x) — Cx. This is a sort of Loren- 
zian to the power of 3/2, which has a width that 
goes to zero when A — > and an amplitude going 
to infmity in this limit. Thus, as A goes to zero, 
the function A°(x) approaches zero for all x except 
for x = rnr/q, n — ...,—1,0,1,.., where it becomes 
singular. Consequently we have 

lim lim A°(x) ^ lim lim A°(x). (46) 

x — >mr A — >0 A — >0 x — >mr 

It is easy to show that the neglected higher order 
term d 2 9 / VA 2 + 4G 2 is independent of the order- 
ing of the limits in eq. I|46|) . Therefor we conclude 
that eq. H4fífl is still vàlid when the neglected term 
is included. In section IVIII we will show that the 
spread of the wave packet causes the adiabaticity 
parameter to be non-zero for other x than x = mr/q 
in the limit A — * 0. 

In fig. [3] we show a plot of one of the màxima 
(x = 0) for A°(x), eq. I)43|l . as function of x and A. 
As expected, the width increase with the detuning 
while the amplitude decrease. It is not clear if sharp 
or broad peaks of the parameter ^4°(a;) favour an 
adiabatic evolution, which will be examined in the 
next section. However, we may note that the integral 
of gSJ 



dxA°(x) 



Po 



mA 



(47) 



is independent of the parameter C, but depends of 
the detuning A. 



VI. NUMERICAL INVESTIGATIONS 

The numerical studies are carried out by simu- 
lating the evolution of wave packets using the split 
operator method [Ï6| . The atòmic wave packet will 
be taken to be in an initial Gaussian bare state in 
its upper level 



$(x,t = 0) = 







(48) 
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FIG. 3: The adiabaticity parameter A°(x) as function 
x and A (scaled units) for the short optical wave length 
situation 11 12> . It is seen that A°(x) becomes singular at 
x — when A — > 0. The other parameters are A = 1, 
q — 1 and po = 2. 



with the Gaussian 



$+(x) 



: exp 



2A? 



exp(ip z)- 



Throughout the paper we use scaled dimcnsionlcss 
parameters in all numerical considerations, such that 
fi = 1 and rn = 1. The 'exact' evolution is obtained 
by evolving (|48|l with the complete Hamiltonian Jïïjl , 
where we will use Case 1 with A = ui, while the adia- 

batic state ^ ad (x,t) 



where 

the states ^^(x,t) are defined in eq. (|30fl . The 
comparison between the two are measured by calcu- 
lating the fidelity 



(50) 



In the case of p « po we choose to plot the fidelity 
as function of (x) ss pot/m + xq, in order to get a 
better idea of where the fidelity changes along the 
mode G(x). 



A. Numerical treatment of a Gaussian mode 

In this section we investigate the adiabaticity cri- 
teria, derived in the previous section, for a Gaussian 
mode shape, while next section studies the periòdic 
mode instead. We have dealt with the limiting cases 
of large and small detunings A, and it was found 
that having an adiabatic evolution supports a small 
coupling, or photon excitation n, in the former, while 
for the latter it is the reverse case. In figs. 01 and 
we display the fidelity HoO|l for moderate and small 
detunings respectively. In the (a) plots the coupling 



strength is larger than in the (b) ones. The pre- 
dicted result is clearly seen; the fidelity is improved 
for small A in fig. 0]and for large A in fig. El In these 
plots we have chosen to plot the fidelity as a function 
of (x) rather than t, which has been done simply by 
transforming (x) ps pqí + xq. This should be justified 
since p is approximately a constant of motion. The 
atòmic wave packet starts at Xq = —200, outside 
the coupling pulse (width is here a = 50), and then 
propagates across the pulse toi» 200 in the first 
example and x = 300 in the second. The reason why 
the packet is propagated longer in the second exam- 
ple is to show how the fidelity goes down when the 
wave packet enters the area, around x = 200, where 
A°(x) has its maximum. This is shown in the inset, 
which gives the location of the maximum (x > 0) as 
function of A for the two coupling strengths. 



(a) 




<*> 



FIG. 4: This and the next figure0show the fidelity l|ïïÜ)l 
as function of (x) and A for the Gaussian mode shape 
® with A = 10 (a) and A = 1 (b). It is seen that the 
adiabaticity/fidelity is increased for smaller couplings A 
when the detuning A becomes large. The other param- 
eters are A^ = 10, xo — —200, po = 5 and a = 50. 



From figs. 0]and[5]we have seen that the expected 
result suggested by the lowest order of approxima- 
tion of the adiabaticity criterion from the previous 
section. However, in the examples p is almost per- 
fectly a constant of motion, which justifies one of 
the approximations. The other one, coming from 
neglecting the higher order term d 2 8, will be consid- 
ered in section |VÏÏ| 



B. Numerical treatment of a standing wave 
mode 

As was pointed out in section for the standing 
wave mode, the limit |A| <C G(x) brought about an 
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FIG. 5: The same as fig. |1J but for small detunings. 
The inset shows the location of the maximum of the 
adiabaticity parameter A(x) as function of the relevant 
detunings and the two coupling amplitudes A = 10 (a) 
and A = 1 (b). The wave packet is propagated to x = 
300 contrary to x — 200 in the previous figure. The 
reason is to see the drop in the fidelity around x ~ 200 
where A(x) has a màxima. The inset indicates that wave 
packets propagating with larger A will 'see' more of the 
maximum of A(x). The other parameters are as in fig. 

H 



adiabaticity parameter A. {%) singular around the 
points x = mr/q with n an integer. At these points 
we can not have |A| -C G(x) as G{x) vanishes there. 

In the optical regime, the atòmic wave packet 
is supposed to extend over several wave lengths, 
and consequently several 5-like non-adiabatic singu- 
lar couplings when A — > 0, and an adiabatic be- 
haviour is not to be expected. The dynamics of wave 
packets evolving according to the Hamiltonian Ipjjl. 
with a coupling g(x) = Asm(qx) and a short opti- 
cal wave length A -C A x has been studied in |13| . 
The periòdic Hamiltonian has a characteristic spec- 
trum of energy bands E v {k) separated by forbidden 
gaps. Here v is the band index and k the quasi mo- 
mentum. In large parameter ranges, including cases 
when |A| <§; A, the atòmic wave packet evolves as 
a 'free' particle with an effective group velocity and 
mass, both determined from the spectrum. Thus, 
it is tempting to expect that the 'freely' evolving 
Gaussian wave packet follows some kind of adia- 
batic behaviour. However, this turns out to be a 
false statement. The 'freely' evolving characteristic 



of the wave packet came from the fact that the initial 
Gaussian bare state l|48|) is approximately an eigen- 
state of the Hamiltonian in these parameter ranges, 
and since the spread of the wave packet is much 
larger than the wave length, the eigenenergy E v (k) 
may be expanded to second order around the initial 
quasi momentum ko. It is a known fact that an ini- 
tial Gaussian wave packet evolving due to an energy 
E = cq + cip + C2P 2 , for some constants co, c\ and 
C2, remains Gaussian. As it turns out, the initial 
Gaussian bare state (|48|) does not approximate an 
eigenstate of the adiabatic Hamiltonian H17fl in the 
small dctuning limit and will therefore not evolve 
as a 'free' particle undcr that Hamiltonain. The fi- 
delity will therefor be small. However, in the large 
detuning limit |A| 3> G(x), adiabaticity is seen, as 
A°(x) <C 1 for all x. The above is shown in fig. |ÜJ 
which displays the fidelity (jf>0p as function of po- 
sition and detuning. For large detunings the adia- 
baticity is seen and the fidelity approaches unity for 
increasing A and decreasing coupling amplitudes A. 




FIG. 6: The fidelity l|50jl as function of (a;) and A for the 
short optical wave length mode shape 142H with A — 1 
(a) and A = 0.1 (b). It is seen that the adiabatic- 
ity /fidelity is increased for smaller couplings A when the 
detuning A becomes large. The other parameters are 
Az = 10 (the width of the wave packet covers several 
wave lengths), Xo = —25, po = 5 and q = 1. 

When À > A x and A — > 0, the atom traverses 
intervals where .4° (2;) <C 1 for the whole extension 
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of the wave packet and the fidelity (|50() should stay 
fairly constant. When the wave packet approaches 
the singular points x = nir/q, a change in the fi- 
delity is expcctcd. For a given detuning and coupling 
amplitude, the width A^ and momentum p roughly 
determine the rate of this change in the fidelity. As 
pointed out in the previous subsection IVI Al it is not 
necessary that a drop in the fidelity is seen when the 
atom traverses a singularity. In fig. [7| we show the 
fidelity for two different coupling strengths A = 1 
(a) and A = 0.1 (b). Here A = 20tt, while = 4. 
It is seen that the main changes of the fidelity oceurs 
at the points x = mr/q. 




FIG. 7: The same as fig. |ü|but for a long wave length; 
A = 20tt and = 4. In (a) A = 1 and (b) A = 0.1. Noto 
that around x = — Wn, 0, IOtt, the fidelity changes the 
most. The initial position and momentum are xo = —50 
and po — 5. 



VII. HIGHER ORDER CORRECTIONS TO 

A t 

In this section we will calculate numerically the 
adiabaticity parameter l|32|l and compare with the 
approximated result of eq. (|35|l . We will also dis- 
cuss the semiclassical time dependent model, which 
is obtained by the replacement p — > po and x — » 
p t/m + x - 

In the approximated result l|35|l . operators where 
replaced by semiclassical parameters and no aver- 
aging over the wave packet was carried out. Thus, 



the atom is seen as a point-like particle moving with 
wcll defined vclocity and position. In order to take 
into account the quantum uncertainties we need to 
average over the wave packet. Note that the aver- 
aging is included already in eq. I|29|l . If averaging 
over the fractions in l|32|) . the relative order of x 
and p must be specified, e.g. (d9)p(A + — A_) _1 ^ 
(A+ - A^)- 1 (d6)p. The separate wave packets 
^^ d {x,t) and ^^ d (x,t) evolves differently, therefore 
At will depend on the initial condition ^(x, 0). The 
effect of averaging depends on the the shapes of the 
wave packets and the coupling, if G(x) various con- 
siderably over the wave packet widths A^ , the aver- 
aging has the effect of smearing out the variations. 
This is the case for a fast oscillating standing wave 
mode in an optical cavity, and therefor this case will 
be left out from this section and we only consider the 
microwave situation for the standing wave mode. We 
also note that if the opposite holds; G(x) is almost 
constant over A x , then we have (G(x)) « G((x)), 
which holds for any smooth function G. 



A. The parameter At for a Gaussian mode 
shape 

The width of the wave packet is assumed to be 
narrower than the width of the Gaussian mode cou- 
pling. This means that the peaks seen in A°(x) 
of fig. Q] should not become consider ably wider 
due to the averaging. The small broadening is ex- 
pected to be more pronounced at later times, that 
is for the second maximum of the adiabaticity pa- 
rameter, since the initially Gaussian wave packets 
"í>^l(x, t) spread with time. In fig. (SJwe give the re- 
sults from numerical calculations of At (sòlid line), 
as defined in eq. (|32")l . together with the approxi- 
mated adiabaticity parameter A°(x) (dashed line), 
where x = pot/m + xq. The circles displays A°(x) 
with the neglected term d 2 9 included. For the up- 
per two plots, the momentum is small (compared 
to A = 100) and therefore varies during the evo- 
lution, while in the lower ones we have p ~ po 
(A = 1), which is also seen by the better agreement 
between dotted/circles and sòlid lines. The left fig- 
ures have a large detuning A = 10 while the others 
have A = 0.05. We note that the main contribution 
to the results comes from the term d9p and not d 2 8, 
as expected. The broadening of the peaks can also 
be seen. 



B. The parameter At for a standing wave 
mode shape 

As already pointed out, when the wave length of 
the standing wave mode is shorter than the width 
of the wave packet the adiabaticity parameter (1:>2I 
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FIG. 8: These figures show the adiabaticity parameter 
1321 1 calculated from wave packet simulations (sòlid line), 
the corresponding approximated result l|35|l (dashed line) 
and the same with the d 2 8 term included in l|35|l (circles). 
In (a) and (b) the coupling amplitude is A — 100 which 
makes p vary considerably, while in (c) and (d) A — 1 
giving that p ~ po . The left plots show the large detun- 
ing situation A = 10, while A = 0.05 for the right ones. 
Note that the dashed and dotted lines are very similar, 
indicating that the main contribution comes from the 
first derivative of 9. The remaining parameters are, (a): 
po = 0.6, x = -200, (b): p = 1.5, x = -300, (c): 
po — 1.5, xo = —200, (d): po = 1.5, xo = —300, and 
a = 50 and A^ = 10. 










2.3 


1 








i 


2 


cx 




>■» 


1 .5 














1 


■1 








Adi 


0.5 











FIG. 9: This shows the same as the previous fig. |H| 
but for the standing wave mode of eq. 1421 . Again, 
the dashed line displays the approximated result l|35|l 
and the circles the same but with the d 2 6 term added. 
The broadening of the peaks, coming from averaging over 
the atòmic wave packet is clearly seen in (b) where the 
peaks are much more narrow. The parameters are, (a): 
A = 0.1, A = 0.5, (b): A = 1, A = 0.05, and for both 
plots A x = 4, xq = —50, po = 5 and q = 0.1. 



will be smeared out due to the averaging over the 
wave packets. Thus, we only consider the case of a 
microwave field, where the extent of the wave packet 
may be smaller than the period of the mode. The 
approximated result (|35|) indicated that for small de- 
tunings, At should be small for all x except at the 
points x — rnr/q where it blowed up. However, the 
averaging will broaden these spikes. This is shown 
in fig. 03 which gives the results of At from wave 
packet propagation, compared with the analytically 
approximated results. The expected broadening is 
especially seen in the plot with narrow peaks, (a) 
A = 0.1, A = 0.5, while in (b) A = 1, A = 0.05. 
In (b) the momentum is not a good constant of mo- 
tion; (p) t ^ po- We note that the main contribution 
comes from the term (d9)p and not d 2 0. 



VIII. RELATION TO SEMICLASSICAL 
TIME DEPENDENT HAMILTONIANS 

In this section we investigate how one may derive 
effective time dependent models by using the adia- 
baticity parameter At- 

Above we loosely argued, that for partides with a 



kinetic energy much larger than the interaction en- 
ergy, we could derive an effective semiclassical time 
dependent model with x = p^t/m + iro- The Hamil- 
tonian then reads 



H = 



I G(t) 



G(t) 



(51) 



where ~ indicates that it is an effective time depen- 
dent model. Here it is understood that t and x are 
related according to x = pot/m + xq. Going to the 
adiabatic frame, by diagonalyzing H with the cor- 
responding unitary operator U as in eq. I|12|) . the 
Schròdinger equation becomes 



"*t(*)" 




'À+ 


+4H 


"%(*)" 


.*!(*). 






.*;(*). 



(52) 

Here the diagonal terms A± are given by the corre- 
sponding expressions H21(l , and the adiabatic correc- 
tion term is 



dt 



o d t e 
-d t è o 



(53) 
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Demanding the off-diagonal term to be small com- 
pared with the 'distance' between the diagonal 
terms, gives the adiabaticity parameter 



Ài 



Sy/ndg(x) fp \ 



(S 2 + 4ng 2 (x)) 



3/2 



m) 



(54) 



x— pot/m+XQ 



which is identical to the approximate result 135|l . 
This has several interesting consequences: 

i) The analysis carried out in this paper concerns 
wave packet propagation, and since the approximate 
adiabaticity criterion we have defined is the same 
as the one obtained from time dependent models, 
much of the conclusions made in this paper are 
applicable to such time dependent models as well. 

it) The constrain At <C 1 is identical with the one 
standardly given in the literature 



(* T (*)lftl*l(*)> 



A+(í)-A_(í) 



< 1. 



(55) 



Recently it has been argued that the criterion 
does not, in general, guarantee adiabaticity 0,] 
The definition IMl'l) takes into account for the initial 
state being used, that is; the atom may start in any 
linear combination of upper and lower atòmic states 
and it is apparent in the population weights N® and 
N?. The corresponding definition for a time de- 
pendent model with instantaneous eigenvalues E n (t) 
and eigenstates \n(t)} would read 



4 = ££i 

n m^n 



(m(t)\d t \n(t)) 



E m {t)-E n {t) 



(56) 



where |c° | 2 is the initial population of state \n(0)). 

in) Deriving an cffective time dependent model as 
above, is only justified when the two wave packets 
^ d (x,t) and ^ d (x,t) follow nearly the same 
trajectory, since x is replaced by a single quantity; 
x = pot/m + xq. It must also be that the mode 
changes smoothly across the wave packet in order to 
neglect d 2 8 and allow (h(x,p)) « h((x), (p)). Both 
these approximations neglect quantum mechanical 
corrections; the term d 2 9 does not show up in the 
semiclassical models since the Schròdinger equation 
is first order in time, but second order in x. The 
quantum uncertainties in x and p are clearly ignored 
in A°{x). 

We now turn to the problem of deriving a time de- 
pendent Hamiltonian when p is not a good constant 
of motion. This method is a kind of inverse problem. 
Given At we may construct a set of Hamiltonians 
H{t) in the form of Q51f) that all give this particular 



adiabaticity parameter. Thus we solve for G(t) from 
the equation 



5d t G{t) 



(5 2 + 4G 2 (í)) 3/2 



= A t , 



with solution 



G(í) = 



/(*) 



(5 2 +4f(t)) 



1/2 • 



(57) 



(58) 



where f(t) = J_ dt' At' ■ This defines a set of time 
dependent Hamiltonians for a given adiabaticity pa- 
rameter, where the coupling G(t) and detuning S 
are related as in l|58|). Note that this set of effec- 
tive Hamiltonians is unique for one particular ini- 
tial condition governing the At- Thus, it gives an 
effective time dependent model for a certain Hamil- 
tonian H{x,p) and initial condition ^>(x,0). For a 
smoothly varying function h(x), over the spread of 
the wave packets, we have (h(x)) w h((x)). Within 
this approximation, and neglecting the d 2 9 term, the 
adiabaticity parameter Ij32|) becomes 



A t = — 



1 

2m 



d6{x í )p i 



A+(x x ) 



ív t ° 



de( X] )p ] 



A+(x T ) 



(59) 

where Xi, Pi, i =[, ] are the expectation vàlues with 
the adiabatic wave functions ^^ Jx, i). When these 
approximations are supposed to hold, it is reason- 
able to expect that we can replace the expectation 
vàlues by their classical counterparts 



dtPi 



dx\ 



d t x\ = 



Pi 



dtPf 



d t xt = 



5A + (x T ) 



dxi 



This is a way to derive a time dependent semi- 
classical model when p is not a constant of mo- 
tion and the wave packets traverse different trajec- 
tòries. The error of such an effective Hamiltonian 
compared to exact ones is left for future research. 
It may be noted that one At defines a set of time 
dependent Hamiltonians of the form l|51|). and it is 
an open question how they are related apart from 
having the same adiabaticity constrains. For ex- 
ample, is there a particular detuning 5 of the time 
dependent model such that the approximated effec- 
tive Hamiltonian becomes an optimal approximated 
model? Note that if the approximation of replacing 
p and x by po and x = pot/m + xo has been carried 
out, we note from l)57p. that we regain the result 
G{t) = G(x)\ x=Pot / m+Xo , provided that we identify 
S = A. As it is defined in the equation, (|51|l is not 
the most general two-level time dependent Hamilto- 
nian, since the diagonal terms are assumed constant. 
If they are allowed to change with time, the method 
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above becomes much more complicated. In general, 
to specify both G(í) and 8 we need, in addition to 
eq. fïïTjl . a second equation. 

IX. CONCLUSION 

In this paper we have discussed the concept of 
adiabaticity within cavity QED models. To do so, 
the original Jaynes-Cummings model, describing the 
atom-cavity field interaction, has been modified by 
taking into account, both the center-of-mass posi- 
tion and momentum of the atom. Two different 
cases have been studied: The atom sees a Gaussian 
or a standing wave mode shape. Numerical simula- 
tions have been used to measure how well the atòmic 
wave packet follows the cavity mode shape adiabat- 
ically. These results have been compared with ap- 
proximated analytically derived criteria for the adia- 
baticity. Good agreement between the two has been 
found. 

The large detvming case, generally called the adia- 
batic limit, implies adiabatic following for any mode 
shape. More interesting is that adiabaticity is, for 
the Gaussian mode, obtained in the opposite limit 
as well, the small detuning case. Contrary to the 
large detuning situation, here a large atom-field cou- 
pling or a high cavity excitation improves adiabatic- 
ity. This is discussed in detail and it turns out that 
for the standing wave mode case, this kind of adia- 



batic limit does not exist. 

The approximations in the derivation of the an- 
alytic adiabaticity criterion have been checked by 
comparing it with the more correctly defined mea- 
sure, which is an average over the atòmic wave 
packet as it traverse the cavity field. The averaging 
tends to smear out ràpid variations of the criterion. 

Finally, we have discussed how the wave packet 
propagation problem may be transformed into an ef- 
fective time dependent semiclassical two-level prob- 
lem. The most simple situation is when p is roughly 
a constant of motion for both of the infernal adia- 
batic wave packets, and one replaces p — ► mv and 
x — ► vt + xo in the original Hamiltonian. The proce- 
dure is not straightforward when p is not a constant 
of motion and the two wave packets traverse differ- 
ent trajectòries. A systematic way to derive a set 
of semiclassical Hamiltonians is presented in these 
cases. It is suggested when such a method is sup- 
posed to work, but no deeper analysis of the validity 
is carried out, and is left for future investigations. 

The paper has been devoted to cavity QED sys- 
tems, but it should be mentioned that the analysis 
can be generalized to a selection of other systems, 
such as: trapped ions, molecular coherent process- 
ing, quantum electrònics, BEC, and also time de- 
pendent models. It is also straightforward to con- 
sider multi-level àtoms and/or multi- modes Jaynes- 
Cummings models. 
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